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It is well known that the hydrodynamics of a zero- temperature superfluid can be formulated in 
field-theoretic terms, relating for example the superfluid four-velocity to the gradient of the phase 
of a Bose-condensed scalar field. At nonzero temperatures, where the phenomenology of a su- 
perfluid is usually described within a two-fluid picture, this relationship is less obvious. For the 
case of a uniform, dissipationless superfluid at small temperatures and weak coupling we discuss 
this relationship within a (p^ model. For instance, we compute the entrainment coefficient, which 
describes the interaction between the superfluid and the normal-fluid components, and the veloc- 
ities of first and second sound in the presence of a superflow. Our study is very general, but can 
also be seen as a step towards understanding the superfluid properties of various phases of dense 
nuclear and quark matter in the interior of compact stars. 
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1. Introduction and summary 

On a microscopic level, superfluidity arises through the formation of a Bose condensate which 
carries charge associated to a ?7 ( 1 ) symmetry of the underlying Lagrangian. The condensed bosons 
may be the fundamental degrees of freedom of the microscopic theory or may be composite objects 
such as Cooper pairs of fermions. In either case, the existence of the condensate allows for friction- 
less transport of the respective charge (mass, in the nonrelativistic context). The reason is that, for 
zero temperature and for sufficiently small velocities of this "superflow", there are no elementary 
excitations that could dissipate energy. This may seem surprising because the Bose condensate 
spontaneously breaks a global symmetry, giving rise to a massless mode - the Goldstone mode - 
that can potentially be excited. However, because the low-energy dispersion relation of the Gold- 
stone mode is typically linear in momentum (and not quadratic or higher order), it turns out that it 
can only be excited beyond a certain critical velocity which equals the slope of the linear part of 
the dispersion. 

At nonzero temperature, on the other hand, excitations of the Goldstone mode (the "phonons") 
are present for any superfluid velocity, and the picture becomes more complicated. Our microscopic 
description of this situation will be as follows. We shall consider a complex scalar field (p with a 
quartic interaction term A<p^. Then we consider the thermodynamics of this system in the presence 
of three independent, external parameters: chemical potential /i, superfluid velocity Vs, and temper- 
ature T . (We will see that pL and are related to the temporal and spatial derivatives of the phase 
of the condensate, respectively. In the two-fluid picture, to be introduced below, our microscopic 
calculation is performed in the rest frame of the normal fluid, i.e., Vs is the velocity of the superfluid 
measured in the normal-fluid rest frame.) We restrict ourselves to a uniform velocity and to small 
temperatures T <^ jj. and small couplings A ^ 1. This allows us to evaluate the effective action, 
the stress-energy tensor etc. analytically. Our approach naturally contains the effect of both the 
condensate and the phonons - and we shall see how both are coupled to each other. It also contains 
in principle the full spectrum of the Goldstone mode, i.e., beyond the linear term mentioned above. 
[For the explicit results we shall truncate the temperature expansion at T^, for which only the linear 
and cubic parts of the dispersion are needed, see Eqs. (3.3).] Due to the nonzero superfluid velocity 
we obtain an anisotropic dispersion relation for the phonons. Although technically a bit involved, 
this approach is completely straightforward from the field-theoretic point of view. The main goal 
is to relate it to the usual "macroscopic" description of a superfluid, that is less common for field 
theorists but widely used in the phenomenology and hydrodynamics of superfluidity. We briefly 
explain this description now. 

On a macroscopic level, a superfluid at nonzero temperatures (below the critical temperature 
where superfluidity breaks down) is usually described in the two- fluid model [1], which was de- 
veloped in the non-relativistic context of superfluid helium. In the two-fluid model, the superfluid 
is viewed as consisting of two fluid components that interpenetrate and entrain each other. Non- 
relativistically, each of the fluid components is characterized by a mass density such that the sum 
of both mass densities adds up to the total (conserved) mass density p = Ps + Pn- While the su- 
perfluid density Ps equals the total density at zero temperature, it decreases with temperature until 
all mass sits in the normal density p,j at and above the critical temperature. One of the predictions 
of this formalism is a second sound mode in superfluids because, besides the usual (first) sound - 
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oscillations in the total density (pressure) - relative oscillations between the two fluid components 
become possible which correspond to oscillations in temperature (entropy). We are interested in 
the relativistic version of this two-fluid formalism [2]. Here, the basic variables ai^e four 4- vectors, 
namely the conserved current j^, the entropy current (also conserved in the absence of dissipa- 
tion), and their conjugate momenta, d^'y^f and @^ (the connection to the microscopic description 
becomes manifest through the phase of the condensate i/a). Two of these 4-vectors are independent, 
i.e., we may for instance choose the two currents as our variables and obtain the momenta as linear 
combinations of them, 

di'\i/ = Bf +As^ , =Af+Cs^. (1.1) 

The coefficients A, B, C depend on the details of the system, i.e., the underlying microscopic theory. 
They can be viewed as derivatives of a generalized energy density A with respect to Lorentz scalars 
built from the cun^ents, see Eq. (4.2). We shall show how to compute these coefficients within our 
scalar field theory. In general, the momenta will not be 4-parallel to their corresponding currents, 
but receive contributions also from the other current. This is interpreted as entrainment between the 
two fluid components and the corresponding coefficient A is called entrainment coefficient. Note 
that here the two fluid components are given by the conserved cun^ent and the entropy current, 
not by the superfluid and the normal fluid. An equivalent formulation [3], somewhat closer to 
the non-relativistic version mentioned above, uses as independent variables one of the currents, 
namely s^, which corresponds to the normal fluid, and one of the momenta, namely d^Y^ which is 
proportional to the superfluid velocity [3]. Both formulations can be translated into each other [4]. 

At first sight, the hydrodynamical variables of the two-fluid formalism seem quite different 
from the variables of the field theory. However, we shall see that there is a one-to-one correspon- 
dence that allows us to compute the basic properties of the superfluid from microscopic physics, 
at least within our simplifying assumptions of a uniform superflow and vanishing dissipation. The 
following points are crucial for this "translation". 

• The field-theoretic calculation is performed in the normal-fluid rest frame. In particular, v.v, 
jj., and T are all measured in this frame. The effective action density can be identified with 
the generalized pressure which in turn is related to the generalized energy density A. 

• With the identification jUv^ = — Vi/a, /i = ^''i//, 7 = 0" and the condition for the normal- 
fluid rest frame s' = 0, eight independent variables are given in the microscopic approach. 
This number equals the number of degrees of freedom of two 4-vectors, for instance the two 
currents j'^ and s^. We shall discuss how the other eight variables s^, &', as well as the 
coefficients A, B, C can be computed as functions of v^, /i, and T. 

In at least two important aspects, this translation between field theory and superfluid hydrody- 
namics is only a first step on which future work should build. Firstly, we only discuss stationary, 
homogeneous cuiTcnts. From the hydrodynamic point of view this is only the simplest situation 
since the hydrodynamic equations d^j^ = d^T^^ = Q with the stress-energy tensor T^^ are triv- 
ially fulfilled in this case. Secondly, since our microscopic calculations are of thermodynamic 
nature, we do not include dissipation. Nevertheless, we may use our results for some nontrivial 
phenomenological features of superfluidity, most notably the two sound modes mentioned above. 
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In particular, we can compute the velocities of first and second sound in the presence of a superflow 
(by superflow we always mean a relative flow between superfluid and normal fluid), see Fig. 1. 

2. Condensate and superfluid velocity 

We start from the Lagrangian for a massless', complex scalar field <p with quartic interactions, 

^ = d^(pd^'(p* -X\(p\'^ , (2.1) 

with coupling constant A > 0. If the chemical potential associated to the conserved U{\) charge is 
nonvanishing, there will be a Bose-Einstein condensate which we parameterize by its modulus and 
its phase, {(p) = -^^^^ ■ A chemical potential can be added by hand, but it can also be introduced 
through the phase of the condensate. In the latter case, the equations of motion allow for the 
following nonzero value of the condensate, 

a 



which suggests that a plays the role of the chemical potential. In the simplest case, usually dis- 
cussed in textbooks, bosons condense in the ground state with vanishing four-momentum. Here we 
allow for a moving condensate, i.e., we want our bosons to condense in a state with non-vanishing 
momentum. This nonvanishing four-momentum is an external parameter and is identical to the 
four-gradient of the phase of the condensate d^\if. In the following, p and are assumed to be 
constant in space-time. By equating the stress-energy tensor from its field-theoretic definition with 
its common hydrodynamic form, one finds that d^Y^^ proportional to the superfluid four- velocity, 

= — ^ . (2.3) 
a 

Therefore, a constant d^'y^f (and thus, via the equations of motion, a constant p) corresponds 
to a static, homogeneous superflow. It is important to understand the meaning of the phase of 
the condensate and of a: with the superfluid three-velocity = — Vi/a/^q'/'^ we can write a = 
doY\/\ — y]. As indicated above, the Lorentz scalar a is the chemical potential, more precisely 
the chemical potential in the frame where the superfluid velocity is zero, which we call the su- 
perfluid rest frame. We see that a differs from d^'y^f only by a factor from a standard Lorentz 
transformation. Hence do^f is the chemical potential in the frame in which the condensate has 
velocity v.,. Below we shall identify this frame with the normal-fluid rest frame in which we de- 
note the chemical potential by /i, i.e., we have pL = doY- We thus see that topological modes of 
the phase play an important role: their winding numbers per unit time and per unit length around 
the U{1) circle determine the chemical potential and the superfluid velocity, respectively. These 
topological modes ai^e different from the topologically trivial oscillations around the ground state, 
which correspond to the collective excitations of the system. By including small oscillations in the 
ansatz for the solution of the equations of motion, one obtains the dispersion relations for the two 
degrees of freedom. In the presence of the condensate, these are the Goldstone mode and a massive 
mode. For small temperatures, only the Goldstone mode becomes populated. We shall now discuss 
the thermodynamics of the system from the microscopic point of view. 

'See Ref. [4] for the case with a nonvanishing mass m. Since many expressions become quite lengthy when in is 
included, we only discuss the case m = in these proceedings. 
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3. Goldstone mode and temperature expansion 



For sufficiently small temperatures (much smaller than the critical temperature) and with the 
assumption of uniformity, the effective action T at the stationary point can be written as^ 

r=-^^-^i:Trln^, S-\k)=( , (3.1) 

with the three-volume V, temperature T, the sum over four-momenta k^ = {ko,k), ko = —i(On with 
the bosonic Matsubara frequencies &)„, and the inverse tree-level propagator in momentum space 
S^^ {k). This expression is obtained by using the zero-temperature result for the condensate (2.2). 
Effects of the melting of the condensate are neglected. The poles of the propagator, i.e., the zeros 
of the determinant of 5^' (^), yield the dispersion relations. Since they are solutions to a nontrivial 
quartic equation, they are very complicated in general. For small momenta, the dispersion of the 
Goldstone mode can be written as 

ek = ci(x)|k| + ^|kp + ... (3.2) 

where we denote x = cos Q with 6 being the angle between the superfluid velocity and the 
momentum k, and where the coefficients in front of the linear and cubic terms are 



v/3-v2(l+ 2x2)yr 3^^2|v,|x 
3-v? 



ci {x) = ^ , (3.3a) 



C2[X) 



(3 - v,^)^ + 2v,^x^(3 - lOvg + 3v^) - (7 - lOv,^ - v^) 
(3 - v?)47rrv?V3-v?(l+2x2) 



3-v?-v^(l+v?) 
(3-v? 



-4|v.|x " . (3.3b) 



The massive mode has a mass £k=o = \/2M \/3 — v? and will not be relevant for the following since 
we assume T <^ jj.. As worked out in Ref. [4], the coefficient of the linear term ci {x) is identical 
to the velocity of first sound, i.e., it can also be obtained from solving the sound wave equations, 
derived from the hydrodynamic conservation equations. The relevant wave equation equation can 
be written as ^^^k^ky = with the so-called "sonic metric" = gf^^ +2v^v^ [6, 7]. More- 
over, one can show that a Lorentz transformation of the four- vector k^^ = (ci (;c)|k|,k) to the frame 
where the superfluid rests yields a transformed velocity c[ = i.e., within our approximation the 
velocity of first sound is isotropic in the superfluid rest frame, even in the presence of a superflow. 

For a translation into hydrodynamics and in particular into the two-fluid framework, we need 
to compute the stress-energy tensor and the cuiTcnt. With the help of the usual field-theoretic 



For temperatures up to the critical temperature a more elaborate, self-consistent formalism is needed [5]. 
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definitions we can write 



1 T 

~X~1V 



-1 



(3.4a) 



(3.4b) 



For tlie stress-energy tensor we have used the gravitational definition, i.e., we take the derivative 
with respect to a general metric and return to g^'^ = (1,-1,— 1,-1) after taking this derivative. 
The evaluation of the effective action (3.1), the current and the stress-energy tensor is now straight- 
forward (up to a renormalization of T^^ , chosen such that the usual thermodynamic quantities in 
the absence of a superflow are obtained [4]). We refer the reader to Ref. [4] for all explicit results. 
As an example, we show the result for the effective action density up to sixth order in temperature. 



y 



AX 



(1 



K^T' (1 



2\2 



10\/3 (l-3v2)2 i05^AtMl-3v2)^ 



.(5 + 30v2+9v^^ 



(3.5) 



For the term, knowledge of the slope of the linear part of the Goldstone mode is sufficient, 
while the term knows about the linear and the cubic parts. Due to the superflow, the stress- 
energy tensor becomes anisotropic, i.e., there is a nonvanishing contribution r"' in the direction 
of V.V, and the diagonal spatial elements have different longitudinal and transverse projections with 
respect to v^, = ^ {8ij — \s,i^s,j) Tij, 7|| = ysj^sjTij. We find that the transverse projection is 
identical to the effective action density (3.5), T± = yF. This is the generalization of the identity 
between the effective action density and the pressure in the isotropic case. In terms of the two-fluid 
picture, the anisotropy of the stress-energy tensor is easily understood: in a single-fluid system one 
can always choose a frame in which the stress-energy tensor is isotropic. For a two-fluid system, 
however, such a frame does not exist in general. In the rest frames of each of the fluid components, 
the pressure in the longitudinal direction with respect to the flow of the other fluid 7|| will differ 
from the pressure in the transverse direction T±. We now explain the two-fluid formalism for a 
nonzero-temperature superfluid. 



4. Covariant two-fluid formalism and sound velocities 

In the superfluid hydrodynamics literature, it is often assumed that the starting point for the un- 
derlying microscopic physics is a "master function" [6] such as the so-called generalized pressure 
*P which is a Lorentz scalar- itself and a functional of the Lorentz scalais a^,&^, and & ■ dy- Each 
of these scalars is in general allowed to depend on the space-time coordinate x (this dependence is 
absent in our simplified field-theoretic treatment above). A Legendre transformation from the con- 
jugate momenta d^yf to the associated currents j^, yields the generalized energy density 
A = — *P + J dY+s ®, and the stress-energy tensor is given by T^^ = —g^^^' + j^d^Y'^^^®^- 
With 

d'V=^, 0'^ = !^, (4.1) 
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one obtains the momenta as linear combinations of the currents (1.1) with the coefficients given by 

dK dA dA 

A = j-—, B = 2j—, C = 2^. (4.2) 
d{j-s) dj^ ds^ 

The connection with field theory is made by realizing that the generalized pressure is identical to 
the effective action density (and therefore identical to the transverse pressure T±), 

*P=^r. (4.3) 

This can be read as a postulate. It also can be derived from the assumption that the field-theoretic 
calculation is performed in the normal-fluid rest frame [4]. This implies that the temperature of our 
statistical ensemble T is measured in the normal-fluid rest frame. We identify T = The relation 
(4.3) requires some explanation since we have just stated that *F only depends on Lorentz scalars, 
while this is obviously not the case for the effective action which depends separately on d^Y^ ^'W^ 
and Since we know how to compute the current [via (3.4a)] and the entropy s = (via its 
thermodynamic definition), we may use Eqs. (1.1) to express A, B, C and the spatial components 
of 0^ in terms of field-theoretic quantities (in the normal-fluid rest frame, the spatial components 
of the entropy current vanish by definition, = 0). The results are 

.-''V^, 5 = -?^, (4.4a, 



where we abbreviated T] = Vv7*'5"i/a+ j " ^^id 



0' ^ 



5V 



.0 , .o.,J-Vv^ 



(VVA)2 



(4.5) 



Again, one can evaluate these expressions in the small-temperature expansion, see Ref. [4] for 
the explicit results. One finds in particular that the entrainment coefficient A vanishes for zero 
temperature which is expected since there is only a single fluid in this case. For any nonzero 
temperature, however, there is, in the language of the two-fluid model, an interaction between the 
two fluid components, characterised by A. One can also translate these coefficients into the normal- 
fluid and superfluid number densities and Hs which are the more natural quantities if the system 
is decomposed in the spirit of the original, nonrelativistic two-fluid formalism, see Eq. (62) of Ref. 
[4]. 

With these relations between field theory and the covariant formalism we can also write down 
the generalized pressure in terms of Lorentz scalars for our particular microscopic model. To this 
end, we compute 



where, in the last two relations, terms of order and are suppressed. This yields, up to order 



t4 ^2 



4A 90V3 



0^ + 2- 



a2 



2 



4 2 
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Figure 1: Polar plots of the velocities of first and second sound ui, for zero temperature T ^ and for 
three different values of the superfluid velocity |Vi | (all velocities measured in the normal-fluid rest frame). 
The direction of the superflow is parallel to the horizontal axis and points to the right; the scale of the axes 
is normalized to the velocity of first sound in the absence of a superflow. 



where we again recover the sonic metric which can be used for a compact notation of the T 
term. In principle, our results also give the term. However, this term is, expressed in terms of 
Lorentz scalars, very complicated and we do not show its explicit form here. 

We conclude with the presentation of the velocities of first and second sound, for which the 
above presented formalism can be used. The complete calculation can be found in Ref. [4]. The 
final results for T — are shown in Fig. 1. Both sound velocities depend on the angle between the 
direction of the sound wave and the direction of the superflow. We find that the velocity of sec- 
ond sound decreases significantly as the superflow approaches its critical value |Vi| = (beyond 
which the system becomes dissipative even at zero temperature). One can show that the velocity of 
second sound receives corrections for small temperatures which increase the velocity. In con- 
trast, there is no temperature correction to the velocity of first sound within our small-temperature 
approximation [4]. 
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